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Abstract. We present the results of a systematic search for multiperiodic pulsators among the Galactic Bulge 
RR Lyrae stars of the OGLE-1 sample. We identify one ’’canonical” double-mode variable (RRd star) pulsating 
in two radial modes. In 38 stars we detect secondary periodicities very close to the primary pulsation frequency. 
This type of multiperiodic variables constitute ~ 23% of RRab and ~ 5% of RRc population of the Bulge. With 
the observed period ratios of 0.95 — 1.02 the secondary periods must correspond to nonradial modes of oscillation. 
Their beating with the primary (radial) pulsation leads to a long-term amplitude and phase modulation, known 
as the Blazhko effect. The Blazhko RRab variables occur more frequently in the Galactic Bulge than in the LMC. 
The opposite tendency is seen in case of the RRd stars. The differences of incidence rates are most likely caused 
by different metallicity of the two populations. We discuss pulsation properties of the OGLE-1 Blazhko stars and 
compare th em w ith predictions of theoretical models. We argue, that the oblique magnetic pulsator model of 
Shibahashi (2000) cannot account for the observations and should be ruled out. 
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1. Introduction 


The Optical Gravitational Tensing Experiment (OGLE; 
Udalski et al. 1992) is devoted to search for dark mat¬ 
ter in our galaxy through the detection of microlensing 
events. As a by-product of this program a large amount 
of photometric data has been collected for stars in the 
Galactic Bulge and Magellanic Glouds. In this series of 
papers we examine short period pulsating variables dis¬ 
covered in the OGLE-1 Galactic Bulge sample (Udalski 
et al. 1994, 1995, 1996, 1997). Paper I (Poretti ^001 ) was 
devoted to monoperiodic S Scuti and RR Lyrae variables. 
In the current paper, second in the series, we will discuss 
multiperiodic RR Lyrae stars. 

The existence of multimode RR Lyrae variables has 
been known since the discovery of double-mode pulsations 
in AQ Leonis (Jerzykiewicz & Wenzel 1977). In follow¬ 
ing years about 300 similar variables have been identi¬ 
fied in several stellar systems, including 181 stars in the 
Large Magellanic Gloud (see Kovacs 2001 for compre¬ 
hensive review). All these variables share common prop¬ 
erties: they pulsate in two radial modes, namely funda- 
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mental and the first overone, with the period ratio of 
Pi/Pq = 0.742 — 0.747. These stars are commonly referred 
to as RRd variables. 

In recent years a new class of multiperiodic RR Lyrae 


variables has been identihed (Olech et al. 1999). They 
are characterised by presence of two or more very closely 
spaced frequencies, with period ratios in the range of 
0.95 — 1.05. Such period ratios are incompatible with the 
radial mode oscillations and point strongly towards pres¬ 
ence of nonradial modes. This new form of multiperiodic¬ 
ity was first discovered in the RRc variables, but later was 
show n to occur in the RRab variables as well (Moskalik 
^000 ). For most c omplete summary of recent observations 
see Kovacs ( 2002 ). 

Massive photometry collected during microlensing sur¬ 
veys is particularly well suited for systematic search of 
multiperiodic RR Lyrae stars. Such a search has already 
been performed in the Large Magellanic Cloud, using 


MACHO data (Alcock et al. 200C; Welch et al. 2002). In 
this paper we present a complete inventory of multiperi¬ 
odic RR Lyrae stars in the OGLE-1 Galactic Bulge sam¬ 
ple. Our primary goal is to establish the incidence rates of 
different types of multimode oscillations in the Galactic 
Bulge population. Comparison between the rates in the 
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Bulge and in the LMC can shed new light on the physical 
conditions favouring such oscillations. We also discuss pul¬ 
sation properties of the identified multiperiodic RR Lyrae 
stars. Preliminary results of this work has already been 
published by Moskalik & Poretti (2002). Here we present 
the complete discussion of our findings. 


2. Search for multiperiodicity 

The OGLE-1 database contains 215 RR Lyrae stars, 150 
of them have been classified as RRab (fundamental-mode 
pulsators) and 65 as RRc (first overtone pulsators). The 
photometric data span ^ 900 days, with typically 130-150 
I-band measurements published for each star. As a first 
step, the lightcurve is fitted with the Fourier sum of the 
form 


mi(t) = (mi)Afccos(27r/cfot-l-<()fe) (1) 

k 

The pulsation frequency fo = 1/Po is also adjusted in the 
fitting process. For several stars slow instrumental drifts 
are present in the data. These are modeled by a cosine 
term with period of 50000 day added to Eq. (1). 

In the next step, the search for additional periodicities 
is performed. This is done with two different, indepen¬ 
dently applied methods: 

1. The fit Eq. (1) is subtracted from the data. The Fourier 
power spectrum of the residuals is then computed in 
order to reveal a secondary frequency, if present. 

2. The fitting formula Eq. (1) is supplemented with an 
additional cosine term with frequency R. The data is 
then fitted for different trial values of R, keeping fixed 
the primary frequency R and the number of its har¬ 
monics, but recalculating their amplitudes and phases 
for each trial. A secondary frequency, if present in the 
data, should reduce the dispersion of the fit in a sig¬ 
nificant way. 

The two procedures give the same results for all the stars 
except one, strengthening our confidence in the frequency 
identifications. In the following, we accept as multiperiodic 
only those variables in which secondary period is detected 
with both methods. 

As a third step, the Fourier fit with two identified fre¬ 
quencies and their linear combinations is performed. To 
this effect, we fit the lightcurve with the following formula: 


mi(t) = (mi) -h ^ Afc„ cos[27r(fcfo -k nRjt -k (j)kn] (2) 

k,n 

The search for additional periodicities is then repeated, 
again using both methods: Fourier transform of residuals 
(method 1) and minimization of dispersion with respect 
to a new trial frequency f 2 (method 2). The process is 
stopped when no significant new terms appear. 

As an example, we show application of the first pro¬ 
cedure described above to the case of BW6 V20. In Fig. |] 


Window 




-.08 -.04 0 .04 .08 


f [c/d] 

Fig. 1. Spectral window for BW6 V20 dataset. 


we present the spectral window of the data. This window 
is typical for OGLE-1 photometry. Fig.^ shows consecu¬ 
tive steps of frequency analysis. Top panel displays Fourier 
power spectrum of the original data. The blow-up of the 
neighborhood of the dominant peak is shown in the right 
column. The second panel displays power spectrum of data 
prewhitened by the first frequency and its harmonics. A 
new frequency has clearly emerged. The next prewhiten¬ 
ing (panel three) reveals another frequency on the oppo¬ 
site side of the dominant peak. Finally, after prewhitening 
by three frequencies no more periodicities are detected 
(panel four). The analysis is stopped. In Fig.|| we display 
the phased lightcurve of BW6 V20, first for the original 
data (top panel), then separately for the three detected 
periods (three bottom panels). 


3. Results 

Most of the RR Lyrae stars in the OGLE-1 sample are 
strictly periodic. Those stars are discussed in Paper 1. 
Departures from a monoperiodic behaviour has been de¬ 
tected in 49 RR Lyrae variables. They can be divided in 
several distinctive types. 
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Fig. 2. Frequency analysis of BW6 V20. Subsequent pan¬ 
els show power spectrum of the original data and the re¬ 
sults of consecutive prewhitening steps. Right column of 
the plot displays fine structure around the main peak. The 
spectra are normalised separately in each panel. 


3.1. RRd stars: variables with radial fundamental and 
first overtone modes 

This type of oscillations has been detected in one object 
of the sample: BW7 V30. The star pulsates with Pq = 
0':‘486889 ± 01*000010 and Pi = 01*36202176 ± 01*00000083. 
The resultant period ratio of Pi/Pq = 0.743541 ±0.000016 
is typical for the RRd variables. The first overtone strongly 
dominates over the fundamental mode, with Aq/Ai = 
0.142. Two harmonics of the first overtone are detected, 
but no harmonics of the fundamental and no frequency 
combinations can be identified. This is perhaps not sur¬ 
prising, considering extremely low amplitude of the fun¬ 
damental mode (Ao = 0.018mag). 


3.2. RR-ivl stars: Blazhko variables with two closely 
spaced frequencies 

In 38 RR Lyrae variables secondary peaks close to the 
dominant pulsation frequency are detected. They are well- 
resolved within our dataset and are not due to a secular 
period variability. Their beating with the primary (radial) 
pulsation results in an apparent long-term amplitude and 
phase modulation, a phenomenon referred to as Blazhko 
effect. 

In most cases (27 stars) only one secondary peak is 
present, forming a close doublet with the primary peak. 
Such frequency pattern occurs in 2 RRc and in 25 RRab 


stars. Following Alcock et al. (200C), we denote these vari¬ 
ables as RRl-:/l and RR0-:/l, respectively. Their main 
properties (period Pq and amplitude Aq of the primary 


BW6 V20 




phase 

Fig. 3. Phased lightcurve (Tmagnitude) of BW6 V20. 
Top: original data phased with the primary period. 
Bottom: lightcurves of the three detected frequencies, fo = 
2.539947c/d, h = 2.554618c/d and fs = 2.525276c/d, re¬ 
spectively. 


peak, frequency separation Af = fi — fg and amplitude 
of the secondary peak Ai) are listed in Table |^. In most 
stars the frequency separation is between 0.008 c/d and 
0.036 c/d, corresponding to a beat (Blazhko) period of 
28 —125 day. The amplitude of the secondary peak is al¬ 
ways very small, never exceeding 0.05 mag. Contrary to 
the case of the LMC stars (Alcock et al. 2000| ; Welch et al. 
p002 ), combination frequencies Mq ± fi = (fc ± l)fo ± Af 
can be detected in about half of the RR-j^l stars of the 
OGLE-1 sample. The amplitudes of all components identi- 











































4 


Moskalik & Poretti: Multiperiodic RR Lyrae variables in the Galactic Bulge 


Table 2. Amplitudes of Fourier components for selected RR-:^1 variables of OGLE-1 sample 


Frequency 

BWC V47 

MM5A V46 

BWll VIO 

BWC V41 

BW3 V13 

BW5 V34 

MM5A V9 

Af 

0.005 

— 

0.014 

— 

— 

— 

— 

fo 

0.067 

0.202 

0.252 

0.221 

0.230 

0.223 

0.251 

fo + Af 

0.050 

0.041 

0.049 

0.041 

0.042 

0.028 

0.016 

2fo 

0.008 

0.068 

0.115 

0.087 

0.103 

0.107 

0.123 

2fo -f Af 

— 

0.028 

0.040 

0.032 

0.027 

0.025 

0.010 

CO 


0.021 

0.072 

0.044 

0.062 

0.072 

0.080 

3fo -f Af 


0.011 

0.033 

0.022 

0.022 

0.024 

— 

4fo 



0.043 

0.028 

0.031 

0.038 

0.052 

4fo -f Af 



0.028 

— 

— 

— 

— 

5fo 



0.027 



0.021 

0.031 

5fo -f Af 



0.024 



— 

— 

6fo 



0.018 



0.016 

0.017 

6fo -f Af 



0.018 



— 

— 

7fo 







0.010 

7fo -f Af 







— 


BWIO V21 

BWIO V20 

BWIO V44 

BWC V61 

BWIO V41 

BWC V51 

BW6 V17 

Af 

fo 

0.136 

0.244 

0.191 

0.096 

0.215 

0.071 

0.146 

fo-f Af 

0.042 

0.023 

0.036 

0.022 

0.049 

0.031 

0.022 

2fo 

0.066 

0.118 

0.083 

0.034 

0.100 

0.019 

0.068 

2fo -f Af 

0.029 

0.022 

0.024 

0.014 

0.032 

0.009 

0.010 

CO 

0.033 

0.062 

0.052 

0.015 

0.052 

0.006 

0.041 

3fo -f Af 

0.023 

0.024 

0.008 

0.016 

0.029 

— 

0.017 

4fo 

0.013 

0.036 

0.028 


0.035 


0.019 

4fo -f Af 

— 

0.016 

0.021 


— 


0.010 

5fo 


0.017 





0.013 

5fo -f Af 


0.015 





— 

6fo 







0.007 

6fo -f Af 







— 


Table 4. Amplitudes of Fourier components for RR-BL variables of OGLE-1 sample 


Frequency 

BW8 V34 

BW6 V20 

BW9 V34 

BWC V15 

BW9 V24 

BW6 V7 

BWl V18 

BWC V33 

BWIO V66 

Af 

— 

— 

— 

— 

— 

— 

0.020 

— 

— 

<1 

1 

A- 

0.030 

0.031 

0.042 

0.072 

0.020 

0.035 

0.032 

0.020 

0.030 

fo 

0.127 

0.244 

0.254 

0.228 

0.206 

0.257 

0.250 

0.175 

0.156 

fo + Af 

0.016 

0.065 

0.063 

0.056 

0.103 

0.042 

0.034 

0.034 

0.053 

2fo - Af 

0.008 

0.028 

0.039 

0.059 

— 

0.034 

0.028 

0.019 

0.027 

2fo 

0.014 

0.104 

0.103 

0.081 

0.068 

0.124 

0.127 

0.087 

0.062 

2fo + Af 

— 

0.042 

0.051 

0.044 

0.066 

0.038 

0.020 

0.033 

0.035 

<1 

1 

A- 

CO 


0.021 

0.030 

0.035 

— 

0.029 

0.019 

— 

0.016 

CO 


0.040 

0.052 

0.036 

0.018 

0.065 

0.080 

0.051 

0.036 

3fo + Af 


0.030 

0.050 

0.024 

0.051 

0.038 

— 

0.030 

0.028 

<1 

1 

A- 


0.012 

0.022 

0.016 

— 

0.024 

0.018 

— 

— 

4fo 


0.022 

0.020 

0.022 

— 

0.038 

0.049 

0.033 

0.017 

4fo + Af 


0.029 

0.027 

0.014 

0.026 

0.030 

— 

0.023 

0.019 

5fo - Af 


— 


0.011 


0.015 

— 


— 

5fo 


0.013 


0.017 


0.018 

0.020 


— 

5fo + Af 


0.015 


0.013 


0.022 

— 


0.014 
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Table 1. RR-z/1 variables of OGLE-1 sample 


Table 3. RR-BL variables of OGLE-1 sample 


OGLE No 

Po 

Af 

Ao 

Al 


[d] 

[c/d] 

[mag] 

[mag] 

RRl-i^l 

BWG V47 

0.256926 

-0.02439 

0.067 

0.050 

MM5A V20“ 

0.391192 

0.04546 

0.126 

0.009 

RRO-i^l 

MM5A V46“ 

0.370187 

0.02530 

0.202 

0.041 

BW3 V17 

0.403331 

0.12005 

0.258 

0.013 

BWll VIO 

0.450658 

0.01864 

0.252 

0.049 

BW3 V13 

0.457517 

0.03255 

0.230 

0.042 

BWG V41 

0.462144 

0.01767 

0.221 

0.041 

BWll V44 

0.482493 

0.03329 

0.204 

0.028 

BW5 V34 

0.490265 

-0.00873 

0.223 

0.028 

BW2 V23 

0.492764 

0.01198 

0.244 

0.044 

BW7 Vis'" 

0.497071 

0.00941 

0.226 

0.035 

BW7 V8 

0.507103 

-0.03312 

0.244 

0.024 

BW7 V33 

0.511226 

-0.00689 

0.261 

0.030 

MM5A V9 

0.516204 

0.03580 

0.251 

0.016 

BWIO V21 

0.521996 

0.01740 

0.136 

0.042 

MM5B V4 

0.524886 

0.02237 

0.184 

0.024 

BWIO V20 

0.527123 

-0.02450 

0.244 

0.023 

BWIO V44 

0.527263 

-0.01206 

0.191 

0.036 

BW8 V15 

0.577234 

0.02543 

0.167 

0.028 

BW5 V36" 

0.594501 

0.01445 

0.139 

0.015 

BW2 V24 

0.597426 

0.01641 

0.118 

0.018 

BWG V61 

0.615946 

0.02145 

0.096 

0.022 

BWIO V41 

0.626882 

0.00855 

0.215 

0.049 

BWl V34 

0.632719 

0.01434 

0.154 

0.031 

BWG V51 

0.649494 

0.01618 

0.071 

0.031 

BW6 V17 

0.651601 

0.00545 

0.146 

0.022 

BWIO V40 

0.682741 

0.01393 

0.117 

0.021 


“ Amplitude and/or phase of primary peak variable. 

^ Second harmonic of primary frequency variable. 
Amplitude and/or phase of secondary peak variable. 


fied in these variables are given in Table & The amplitudes 
and frequencies listed in Tables ^ and g are determined 
through the least-squares fit of Eq. (2) to the data. 


3.3. RR-BL stars: Blazhko variables with three 
symmetrically spaced close frequencies 

In several RR Lyrae variables two secondary frequencies, 
located on opposite sides of the primary peak, are present. 
Together with the primary frequency they form a close 
equally spaced triplet. Such frequency pattern is found 
in 1 RRc and in 9 RRab stars. Eollowing Alcock et al. 
( 200C ), we denote these variables as RRl-BL and RRO- 
BL, respectively. Their properties are listed in Table 
The amplitudes of secondary peaks are somewhat larger 
than in the RR-z/1 stars, reaching up to 0.1 mag. The fre¬ 
quency separations Af are similar to those observed in 
RR-z/1 stars. With the exception of one star, the triplets 
are exactly equidistant. This has been verified by fitting 
the data with Eq. (2), with all 3 frequencies of a triplet 


OGLE No 

Po 

[d] 

Af_ 

Af+ 

[c/d] 

Ao 

[mag] 

A_ 

A+ 

[mag] 

RRl-BL 

BW8 V34 

0.320365 0.01544 

0.127 

0.030 





0.016 

RRO-BL 


BW6 V20 

0.393709 

0.01467 

0.244 

0.031 

0.065 

BW9 V34 

0.449297 

0.00896 

0.254 

0.042 

0.063 

BWG V15 

0.458717 

0.00376 

0.228 

0.072 

0.056 

BW9 V24 

0.476339 

-0.00185 

0.00239 

0.206 

0.021 

0.103 

BW6 V7 

0.525015 

0.03363 

0.257 

0.035 

0.042 

BWl V18 

0.529560 

0.01876 

0.250 

0.032 

0.034 

BWG V33 

0.550304 

0.00865 

0.175 

0.020 

0.034 

BW6 V29 

0.562906 

0.02283 

0.164 

0.019 

0.022 

BWIO V66 

0.581234 

0.01410 

0.156 

0.030 

0.053 


Table 5. RR0-:^2 variable of OGLE-1 sample 


OGLE No 

Po 

Af 

Ao 

Al 


[d] 

[c/d] 

[mag] 

[mag] 

BWG V106 

0.464967 

0.00708 

0.02675 

0.188 

0.020 

0.020 


treated as independent parameters. The resulting sepa¬ 
rations —Af_ = — (f_ — fo) and Af+ = f+ — fo dif¬ 
fer by less than 0.00005c/d, which within our dataset 
is indistinguishable from zero. Thus, in the final solu¬ 
tion Af_ = — Af+ is assumed. The one deviating case 
is BW9 V24, where Af_ -b Af+ = 0.00054 c/d. We note, 
that in this particular star the triplet is barely resolved. 
Since the A_ component in BW9 V24 is very weak, the 
insufficient resolution can have a significant effect on its 
frequency determination. The deviation from equidistant 
spacing seen in BW9 V24 needs to be confirmed with a 
longer dataset. 

In all but one RR-BL stars we detect combination fre¬ 
quencies fcfo ± Af. The amplitudes of all identified com¬ 
ponents are given in Table ^ The combination peaks, to¬ 
gether with harmonics of fo, form a sequence of equidistant 
triplets. This structure can be followed up to the third or 
fourth harmonic. We note in passing, that the observed 
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Fig. 4. Seasonal lightcurves (I-magnitude) of suspected Blazhko variables. Panels from top to bottom correspond to 
observing seasons of 1992, 1993 and 1994, respectively. 


Table 6. Variable period RRc stars of OGLE-1 sample 


Table 7. Suspected Blazhko variables of OGLE-1 sample 


OGLE No 

Po 

[d] 

Ao 

[mag] 

P 

[s/s] 

P/p 

[yr] 

BWG V37 

0.380122 

0.140 

1.7x10"'^ 

6.1x10® 

BW2 VIO 

0.507782 

0.118 

3.1x10"'^ 

4.5x10® 

BW7 V51 

0.272148 

0.154 

2.9x10"® 

2.5x10'* 

BW9 V38 

0.305731 

0.117 

-4.0x10"® 

-2.1x10* 

BWll V55 

0.286179 

0.123 

-5.1x10"® 

-1.5x10* 

MM5A V24 

0.379396 

0.093 

-5.4x10"® 

-1.9x10* 


frequency pattern is the same as found in several well- 
studied field RR Lyrae B lazhko v ariabl es (Bor kowsk i |l980| ; 
Smith et al. |l994 , 199E ; Kovacs 1995[ Nagy |1998| Szeidl 
& Kollath 120001) . 


3.4. RR-u2 stars: Blazhko variables with three 
nonequidistant close frequencies 


OGLE No 

Po 

Al 

4>i 


[d] 

[mag] 

[rad] 

RRc 

BWll V34 

0.341674 

0.139 ±0.002 

0.00 ± 0.02 



0.122 ±0.002 

0.14 ±0.02 



0.119 ±0.003 

0.03 ± 0.02 

RRab 

BW4 V5 

0.474673 

0.167 ±0.013 

0.00 ± 0.05 



0.223 ± 0.004 

0.13 ±0.02 



0.182 ±0.003 

-0.04 ± 0.02 

BW5 V13 

0.494952 

0.065 ± 0.006 

0.00 ±0.10 



0.326 ± 0.006 

-1.42 ±0.04 



0.286 ± 0.005 

-1.01 ±0.02 

BW6 V18 

0.541402 

0.290 ± 0.014 

0.00 ± 0.03 



0.199 ±0.004 

-0.22 ± 0.02 



0.205 ± 0.005 

-0.24 ± 0.03 


Yet another pulsation patter is found in BWC V106. In 
this star, two secondary frequencies are detected, but both 
located on the same side of the primary peak. Together 
with the primary frequency they form a close triplet, 
which is neither equally spaced, nor centered on the pri¬ 
mary. Following Alcock et al. ( pOOC ) we clasify BWG V106 
as RR0-^2-type variable. This type of pulsation pattern 
is extremely rare, only 3 similar stars have been identi¬ 
fied among 1327 first overtone RR Lyrae pulsators in the 
LMG (Alcock et al. 2000| ). The properties of BWG V106 


are listed in Table ||. No combination peaks are detected 
in this star. 

3.5. Miscellaneous RR Lyrae variables 

In 10 RR Lyrae stars we find after prewhitening a sig¬ 
nificant residual power at frequency almost identical to 
that of the (just removed) primary component. Such be- 
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Table 8. Variable Types in OGLE-1 Galactic Bulge Sample: First overtone RR Lyrae Stars (RRc Stars) 


Type Description Number % LMC % 

RRc . single-mode 54 83.1 ± 11.3 69.0 ±1.3 

RRd . F/IH double mode 1 1.5 ± 1.5 13.6 ± 0.9 

RRl-i^l . two close components 2 3.1 ± 2.2 1.8 ± 0.4 

RRl-BL . symmetric triplet 1 1.5 ± 1.5 2.1 ±0.4 

RRc-PC . period change 6 9.2 ± 3.8 10.6 ± 0.9 

RRc-NC. nonclassified 1 1.5 ± 1.5 0.5 ± 0.2 


Table 9. Variable Types in OGLE-1 Galactic Bulge Sample: Fundamental Mode RR Lyrae Stars (RRab Stars) 


Type Description Number % LMC % 


RRab . single-mode 112 74.7 ± 7.1 ? 

RRO-i^l . two close components 25 16.7 ± 3.3 6.5 ± 0.7 

RRO-J/2 . three close components 1 0.7 ± 0.7 ? 

RRO-BL . symmetric triplet 9 6.0 ± 2.0 3.7 ± 0.5 

RRab-PC . period change - - ? 

RRab-NC . nonclassified 3 2.0 ± 1.2 ? 


haviour is a signature of a slow phase and/or amplitude 
variability, not resolved within the lenght of available data 
(~ 900 days for OGLE-1). 

For closer examination of the long-term variations 
in these objects we have performed Fourier fits of their 
lightcurves separately for each observing season. In 6 stars, 
all pulsating in the first overtone, we find the pulsation 
phase to vary from year to year, but the amplitude of 
the pulsation remains constant. These stars are listed in 
Table |[ Glearly, they undergo a secular period change. 
Assuming, that the period varies linearly in time, we 
can estimate for each object the slowest dP/dt consistent 
with the data. We show them in the last two columns 
of Table ^ We stress, that the numbers derived here are 
based on three observing seasons only and therefore should 
be treated as preliminary. Nevertheless, we can conclude 
that in all 6 cases pulsation period changes on a timescale 
of no more than ^ 10^ yr, which is much shorter than the 
timescale predicted by the theory of stellar evolution. 

In the remaining 4 stars (1 RRc and 3 RRab-type) both 
amplitude and phase vary from year to year. We list these 
objects in Table |^, where in the last two columns we give 
the values of the Fourier amplitude Ai and phase (see 
Eq. (1)) for each observing season. Seasonal lightcurves of 
the stars are displayed in Fig.^. The amplitude variability 
is hardly noticable in the RRc star BWll V34, but in the 
RRab stars it is rather obvious. In the extreme case of 
BW5 V13 the amplitude increases by factor of five within 
one year. The change of the pulsation amplitude is ac¬ 
companied by change of phase and change of shape of the 
lightcurve . This last effect is most noticable in BW4 V5 
and BW6 V18. All 4 stars are most likely Blazhko variables 
with modulation period (beat period between the modes) 
longer than the span of available data. However, with the 
data in hand we are not able to make any firm statement 
about the nature of the observed long term behaviour. 


In the following discussion we will treat all 4 variables as 
” nonclassified”. 


4. Discussion 

4.1. Incidence rates 


Tables g and g present the inventory of different types of 
RR Lyrae variables in the OGLE-1 Galactic Bulge sam¬ 
ple, separately for the overtone and for the fundamen¬ 
tal mode pulsators. The canonical double mode variables 
(RRd stars) are somewhat arbitrarily included into the ta¬ 
ble for overtone variables. This is justified by strong dom¬ 
inance of the overtone mode in these stars - as a result 
they are usually disguised, until full frequency analysis, as 
RRc stars. For each type of variability we list the number 
of stars found (third column ) and the estimated incidence 
rate in the population (fourth column). The standard de¬ 
viations of the incidence rates are calculated assuming 
Poisson distribution. For the purpose of comparison, we 
also show the corresponding rates for the LMG sample 


(Alcock et al. 2000; Welch et al. 2002), these are listed in 
the last column of the tables. The LMC rates have not 
been published for all the subclasses; the missing entries 
in the tables are substituted with the question marks. 

Several important conclusions can be drawn from data 
of Tables and 0: 


1. Blazhko variables (RR-i/1, RR-z/2, RR-BL) occur 
more frequently among fundamental-mode pulsators 
(RRab) than among overtone pulsators (RRc). The 
property is common for both Galactic Bulge and LMC 
populations, being somewhat stronger in the Bulge. 

2. For the RRc stars, the fraction of Blazhko variables 
is (within statistical error) the same in the Galactic 
Bulge and in the LMC. In case of the RRab stars, 
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however, the fraction is twice higher in the Bulge 
(23.3 ± 3.9% vs. 10.2 ± 0.8%). 

This difference does not result from lower noise level 
of the Galactic Bulge data. The weakest secondary 
peaks detected in the LMC stars have amplitudes of 
~ 0.020 mag. Adopting the same threshold for our 
sample, we are still left with 31 RRab Blazhko stars, 
which is 20.7 ± 3.7% of the RRab population. We con¬ 
clude, that the higher incidence rate of RRab Blazhko 
stars in the Galactic Bulge, as compared to the LMG, 
is real. It is tempting to speculate that this effect is 
caused by difference in metallicity of the two popula¬ 
tions. 

3. The RRd variables are significantly less common in the 
Galactic Bulge than in the LMG. Assuming the same 
incidence rate as in the LMC, we would expect to find 
at least 7 variables of this type among 65 RRc stars of 
the Bulge (3cr lower limit). Instead, we find only one. 
Again, different metallicity of the two populations is 
the likely culprit. 

4. Secular period change is detected in about 10% of all 
Galactic Bulge RRc stars. This is the same fraction as 
for the RRc stars in the LMC. Interestingly enough, no 
RRab stars in the Bulge show this type of behaviour. 
We remind here, that the observed period variability 
is too fast to be of evolutionary nature. Its origin is 
currently unknown. 


4.2. Pulsation properties of RRab Blazhko stars 

With the sample of 35 fundamental mode Blazhko stars, 
we are in position to discuss the group properties of 
this type of variables. In Fig.^ we display the Period- 
Amplitude diagram for the Galactic Bulge RRab stars. 
In the upper panel we plot peak-to-peak amplitude, as 
given by OGLE-1 catalog. In the lower panel we plot the 
Fourier amplitude of the radial mode, as determined by 
the Fourier fit Eq. (1) or Eq. (2). Monoperiodic variables 
are marked by asterisks, the Blazhko variables by open cir¬ 
cles (RRO-zzl and RRO-i/2) and by filled circles (RRO-BL). 
It is immediately apparent that the presence of secondary 
frequencies neither depends on nor affects the amplitude 
of the primary (radial) pulsation. It does not depend on 
the pulsation period, either. The stars with close frequency 
doublets (RRO-z^l variables) are detected with roughly the 
same probability at all periods represented in the sample. 
The occurrence of RRO-BL variables, however, seems to 
be limited to P < 0‘?6. Nevertheless, the overall fraction of 
Blazhko variables is essentially period independent. This 
conclusion is also supported by the histogram shown in 
Fig.|. The fraction of Blazhko variables in each period bin 
is within Icr consistent with the average value of 23.3%. 

In Fig. ^ we show the frequency separation Af = f^ — fp 
for RRab Blazhko variables of our sample. Af displays no 
trend with the radial mode period Pq. For 80% of RRO- 
vl stars frequency separation is positive, corresponding to 
secondary frequency being higher than the primary one. 



Fig. 5. Period-Amplitude diagram for OGLE-1 RRab 
stars. Top: peak-to-peak amplitudes from OGLE-1 cat¬ 
alog, bottom: Fourier amplitudes of the radial mode. 
Monoperiodic variables marked by asterisks, Blazhko vari¬ 
ables by open circles (RRO-J^l and RR0-zz2) and filled cir¬ 
cles (RRO-BL). 
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Fig. 6. Period distribution of OGLE-1 RRab stars. Black 
area corresponds to Blazhko variables. 



Identical distribution of Af has also been found for RRab 
stars in the LMG (Welch et al. 2002). The negative values 
of Af occur only in the narrow range of periods between 
0'?49 and 0'?53. Interestingly, this particular period range 
is avoided by RRO-BL stars. The values of Af for triplets 
(RRO-BL stars) are on average slightly smaller than for 
the doublets (RRO-izl stars). In 8 out of 9 RRO-BL stars 
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Fig. 7. Frequency difference, Af, for OGLE-1 RRab 
Blazhko variables. Same symbols as in Fig.||. 


the higher amplitude secondary peak has also higher fre¬ 
quency. Again, very similar statistics has been found for 
the RRab stars in the LMC. 

While Af distribution for the RRab stars is very sim¬ 
ilar in both galaxies, it is distinctively different than for 
the RRc stars of the LMC (Alcock et al. 200C ). In the 
latter case, the frequency separation is negative in more 
than half of the doublets (RRl-z/1 stars) and the higher 
amplitude secondary peak in the triplets (RRl-BL stars) 
can occur on either side of the primary peak, with about 
equal probability. The values of |Af| in the LMC RRc 
stars are also noticably larger than in the RRab stars. 
Unfortunately, the number of RRc Blazhko variables in 
the Galactic Bulge is too small for a meaningfull discus¬ 
sion. Nevertheless, even for this very limited sample we 
find the same result: average Af is larger than for the 
RRab stars. Apparently, whatever mechanism is respon¬ 
sible for excitation of the secondary frequencies, it works 
diferently in the overtone and in the fundamental mode 
RR Lyrae pulsators. 


4.3. Blazhko variables and metallicity 

In Fig.^ we display the Fourier phases (/) 3 i and ^41 for all 
the Galactic Bulge RRab stars. This plot for the monope- 
riodic pulsators has already been presented in Paper I 
(their Fig. 12), here we have added the multiperiodic 
(Blazhko) stars. In case of multiperiodic stars, the Fourier 
phases are derived for the primary (radial) component of 
the pulsation. Specifically, they are constructed with the 
phases of the fcfo components of the Fourier fit Eq. (2). 
We stress that this is a different procedure than used by 
Jurcsik et al. ( 2002| ), who have evaluated the Fourier pa¬ 
rameters of the snapshot lightcurves at various phases of 
the beat (Blazhko) period. As a result, their parameters 
are time dependent. The phases derived in the current 
paper do not vary in time. They describe the shape of 
the primary (radial) pulsation after secondary frequencies 
have been filtered out. 



Fig. 8 . Fourier phases (fsi and (/) 4 i for OGLE-1 RRab 
stars. Same symbols as in Fig. |[ 


In the case of monoperiodic stars the Fourier phases 
(/) 3 i and (/) 4 i follow a single progression up to P ^ 0'?55. At 
longer periods the progression splits into three separate 
tails. As has been shown in Paper I, these tails are formed 
by stars of different metallicity, with the upper tail corre¬ 
sponding to the highest and the lower tail to the lowest 
value oHFe/H]. 

Fig.N shows that the Blazhko variables of our sample 
choose preferentially the upper tail of the progression. This 
is best visible in the plot of ((> 31 . Out of 11 Blazhko stars 
with P > 01*55, 7 fall on the upper tail, while only 2 belong 
to the central and to the lower tails, each. In the ^41 plot, 
these variables are split between the tails in the same way. 
The distribution of Fourier phases implies that the inci¬ 
dence rate of Blazhko variables increases with metallicity. 
We note, that the same dependence has already been sug¬ 
gested by comparison between the Galactic Bulge and the 
LMC (see Sect. 4.1). 

5. Conclusions 

We have conducted a systematic search for multiperi¬ 
odic pulsators in the OGLE-1 sample of Galactic Bulge 
RR Lyrae stars. Multiperiodicity has been established for 
39 variables. In addition, we have identified 6 RRc stars 
with changing pulsation period and 4 stars exhibiting 
long-term change of pulsation amplitude. 

Among 39 multiperiodic variables only one is a “canon¬ 
ical” double-mode pulsator (RRd stars), with two radial 
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inodes - fundamental and first overtone - excited. In the 
remaining 38 variables a different type of multiperiodicity 
is found - additional peaks are present very close to the 
primary pulsation frequency. Their beating with the pri¬ 
mary (radial) pulsation results in an apparent long term 
amplitude and phase modulation. This phenomenon is re¬ 
ferred to as Blazhko effect. Majority of the Blazhko vari¬ 
ables come in one of two flavours: we detect either a single 
secondary peak, forming a doublet with the primary fre¬ 
quency (RR-z/l stars), or a pair of secondary peaks, which 
together with the primary frequency form an equidistant 
triplet centered on the primary peak (RR-BL stars). While 
the frequency triplet can result from periodic amplitude 
and/or phase modulation of a purely radial pulsation, such 
a process cannot produce a doublet. The observed period 
ratios of Pi/Po = 0.95 — 1.02 are not compatible with the 
excitation of two radial modes. This implies that a sec¬ 
ondary component of the doublet must correspond to a 
nonradial mode of oscillation. 

The Blazhko variables occur among both RRc (over¬ 
tone) and RRab (fundamental-mode) stars, but the inci¬ 
dence rate is not the same. In case of the Galactic Bulge 
it is 5% for the RRc stars and 23% for the RRab stars. 

The incidence rates of multiperiodic RR Lyrae vari¬ 
ables differ between the Galactic Bulge and the Large 
Magellanic Gloud. Specifically, the double-mode (RRd) 
stars occur much more frequently in the LMG, while the 
opposite is true for the RRab Blazhko stars. We believe 
that the effect is caused by different metallicity of the 
two populations. Such interpretation is supported by the 
Fourier analysis of the Bulge RRab stars. The hypothesis 
of metallicity influence can be further tested when statis¬ 
tics of multiperiodic variables in the SMG is established. 

The occurence of Blazhko variability among RRab 
stars correlates neither with amplitude nor with period of 
the primary pulsation. In vast majority of doublets (RRO- 
vl stars) the secondary frequency is higher than the pri¬ 
mary one. The same is also true for the higher of the two 
secondary peaks in the triplets (RRO-BL stars). This is 
a different distribu tion than observed in case of the RRc 
stars (Alcock et al. 200C ). The values of the frequency sep¬ 
aration Af are also different, being systematically smaller 
than in the RRc stars. Clearly, the detailed picture of the 
Blazhko variability depends on the primary mode of pul¬ 
sations. 

The observational properties discussed in the current 
paper put constraints on any proposed explanation of 
the Blazhko effect. Two models are most popular nowa¬ 
day: the oblique magnetic pulsator model (Shibahashi 
1995 , 2000| ) a nd th e 1 1 resonance model (Nowakowski & 


Dziembowski 2001). The results presented here seem to 
rule out the first of the two. The oblique pulsator model 
predicts splitting pulsation frequency of a mode into an 
equally spaced quintuplet. We never detect such a struc¬ 
ture, despite specifically looking for it. We see only triplets 
or doublets. The model cannot explain a striking differ¬ 
ence of the Blazhko incidence rate between RRc and RRab 
stars. If the modulation is indeed due to presence of the 


magnetic field, then it should occur equally likely inde¬ 
pendently of the pulsation mode. Finally, as the pulsation 
amplitude is supposed to vary with the rotation period of 
the stars, it is hard to understand why Blazhko periods 
of RRc stars are systematically shorter than that of the 
RRab stars. 

The resonant model avoids all of the above difficul¬ 
ties. It generates equally split triplets, not quintuplets. 
It can explain naturally the asymmetry of amplitudes in 
the triplet. The model predicts higher incidence rate of 
Blazhko variables amon g RR ab than among RRc stars 
(Dziembowski & Cassisi 1999 ). This is in qualitative (but 
not quantitative !) agreement with the observations. The 
modulation period is not directly related to the rotation 
period of the stars, but is determined by interplay between 
the frequency spacing and rotational splitting of the non¬ 
radial modes involved in the interaction. As such, it can be 
different for the overtone and for the fundamental-mode 
pulsators. 

The outstanding difficulty of the resonant model is 
explaining the existence of frequency doublets. In some 
cases, the observed doublet can actually be a triplet, with 
one sidepeak too small to be detected. For several stars, 
however, such an explanation would require an extreme 
amplitude asymmetry in order to be valid. For example, 
we find 3a upper limits of A+ /A_ < 0.08 for BWG V47 
(RRc star) and A_/A+ < 0.19 for BWG V51 (RRab 
star). It is not clear at present if so strongly asymmet¬ 
ric triplets (or pure doublets) can be reproduced by the 
resonant mode coupling theory. 
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